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Abstract: We look for ways to destabilise the vacuum. We describe how dense matter en- 
vironments source a contribution to moduli potentials and analyse the conditions required 
to initiate either decompactification or a local shift in moduli vevs. We consider astro- 
physical objects such as neutron stars as well as cosmological and black hole singularities. 
Regrettably neutron stars cannot destabilise realistic Planck coupled moduli, which would 
require objects many orders of magnitude denser. However gravitational collapse, either in 
matter-dominated universes or in black hole formation, inevitably leads to a destabilisation 
of the compact volume causing a super-inflationary expansion of the extra dimensions. 
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1. Introduction 

String theory has no free parameters. All coupling constants are instead determined as 
vacuum expectation values of scalar fields - moduli. The values these scalar fields take 
are determined by the moduli potential, and these values determine the parameters of the 
Standard Model and through them the masses, couplings and interactions of all known 
particles. 

Moduli potentials have many ingredients, and much work has been done on construct- 
ing potentials that stabilise the moduli in phenomenologically attractive fashions. However 
moduli vevs are environmental and so there is no reason in principle why they should sta- 
bilise at the same values at all regions in space and time - indeed, we should expect the 
converse. In this article we therefore look at ways of destabilising moduli from their vacuum 
values. 

The basic mechanism we investigate is simple: as moduli source the couplings of Stan- 
dard Model fields, any form of energy density represents a source for the moduli potential. 
If the local energy density is sufficiently large - and large energy densities are realised both 
within neutron stars and in the context of cosmological singularities - then this gives a con- 
tribution to the moduli potential that can destabilise the modulus vev from its minimum. 

Destabilised moduli vevs may be hard to achieve but the payoff if it can be done is 
large: 
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• Small changes in moduli vevs would give a continuous deformation away from the 
Standard Model, with associated small shifts in particle masses and couplings. 

• Any region in which moduli vevs are altered can catalyse exotic processes that are 
forbidden or highly suppressed within the Standard Model. A simple example is that 
of proton decay. Within the Standard Model the proton can decay via electroweak 
sphaleron processes. However these are nonperturbative and are suppressed by A ~ 

fS7T 

exp(— S) ~ e 9 sc/(2) anc i tn e resulting proton lifetime is of order 10^ years with 
J\f ~ 100. Any local region in which the SU(2) coupling were substantially stronger 
could therefore catalyse the decay of protons that enter it. 

• Large shifts in moduli can potentially provide windows into entirely different vacua 
of the underlying theory, which could perhaps be realisable locally as stable solitonic 
objects. 

For these reasons we think it worthwhile to investigate the possibility of destabilising 
moduli in various contexts. This discussion can be made concrete by recent developments in 
moduli stabilisation, as it is difficult to discuss the chances of destabilising moduli without 
concrete and well-motivated potentials that first stabilise moduli. 

For other related work, see for example [1, 2, 3, 4]. 



2. The effective field theory 

We study the possibility and the consequences of moduli/matter interaction within the 
framework of four dimensional supergravity. Starting from the Type IIB action in 10 
dimension, upon dimensional reduction in the large volume limit one finds that the action 
for the moduli sector is 



S = f d^x^-g (Gfjd^fd^ - V(<f>, 0)) , (2.1) 



where Gq = didjK is the Kahler metric and (pi are the moduli fields of the theory, i.e. the 
dilaton, complex structure and Kahler moduli. The F-term potential takes its usual form 

V(<f>, 4>) = e K (G fj D i WD ] W - 3\W\ 2 ), (2.2) 

where DiK = diK + KdiW. 

We will work within the LARGE volume models [5], which allow for the stabilisation 
of all the moduli while generating a nonsupersymmetric AdS minimum at exponentially 
large volumes. We briefly review the properties of these models and the masses of the 
moduli in them. These models are characterised by the inclusion of o/ 3 corrections to the 
Kahler potential [6]: 

/ £ (S±S)3/2\ 

K = -2\n [V+ ^ V \+K cs +K s , (2.3) 
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as well as nonperturbative effects in the superpotential W [7]: 

W = W Q + s £ j A i e ia ^. (2.4) 

i 

It was shown in [5] that these models have AdS minima with broken SUSY for ge- 
ometries in which one Kahler modulus controls the overall volume while the remaining 
Kahler moduli describe collapsible four cycles. This so called 'Swiss cheese' geometry is 
quite generic and various explicit cases have been described in [8, 9, 10]. In this work we 
will focus on the simplest realisation of this geometry, the two modulus model defined by 
a degree 18 hypersurface in the weighted projective space Pj^ 1 1 6 g j . The volume for this 
orientifold is given by 

where the modulus t\, parametrizes the volume of the extradimensional manifold and t s 
the size of a 'small' four cycle in the manifold. We note that the results found here should 
remain qualitatively the same for more complicated manifolds provided they are of the 
Swiss cheese type. 

The scalar potential, Eq. (2.2), for the LARGE volume models can be computed by 
taking the Kahler potential and the superpotential, Eqs. (2.3) and (2.4). One can choose 
to stabilise the axio-dilaton and the complex structure moduli in a SUSY preserving way, 
DW = , by turning on fluxes in the extra dimensions. This will cause their contribution 
to the scalar potential to vanish. For the two moduli model we are considering the scalar 
potential can be written as: 

-2ar s 



where 



V = C^ C 2 T -^_ + % (2.6) 



& = -Xa 2 \A\ 2 , (2.7) 

C 2 = 4\AW\a, (2.8) 

3 \W\ 2 £ 

4 9s' 



This potential has a LARGE volume AdS minimum at 



2gl' 2 V 2ar s 

3\W\Jfse aTs ( 3 
V= ' ,.| (l-i • (2-11) 



Aa|^4| V 4ar s 
2.1 Moduli spectrum 

Given that the moduli have been stabilised at the minimum of the potential one can 
compute their masses: 

™ 2 _ n fi 

dr? 



-r) 2 V 

m 2 = G^—^ (2.12) 
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where in the particular case under consideration i £ {s, b}. Noting that the leading terms 
in the volume expansion of the inverse Kahler metric are 



G ss " oc V, (2.13) 
G bl oc V 4/3 , (2.14) 

once we compute the second derivatives of the scalar potential Eq. (2.6) and restore the 
factors of Mp we find 

m s «y, (2.15) 

m b «^. (2.16) 

In the LARGE volume minimum one finds that the small modulus is much heavier 
that the large modulus. One can therefore consider that it decouples from the theory, and 
take T;, to be the only dynamical variable in the problem. This corresponds to using 

e ar s = y(i + su bleading terms), (2.17) 

to eliminate t s dependence. We will later apply this to the scalar potential, Eq. (2.6), but 
first let us formulate the problem in terms of canonically normalised fields. 

2.2 Canonical normalisation of the volume modulus 

As discussed in the previous section, the moduli associated with the volume of the 'small' 
four cycle gets, through moduli stabilisation, a very large mass and can therefore be in- 
tegrated out. What remains is the theory of a single scalar field. The kinetic part of the 
Lagrangian is 1 

L K = GadpTtdPTi = -^d^n. (2.18) 
For convenience we will work with the canonically normalised field defined by 

±^ d ^ Tb= l - m \ (2.19) 

so we find 

$ = y|lnr 6 = y|lnV. (2.20) 
Using Eqs. (2.17) and (2.20) one rewrites the potential, Eq. (2.6), as 

27<J> 



V = (l-a$ 3/2 )e~V^ <I> , (2.21) 

where we have ignored factors of order unity (note this requires including the subleading 
terms of eq. (2.17)). In Fig. 1 we plot the potential for the canonically normalised volume 
modulus, Eq. (2.21), which exhibits the characteristic AdS minimum at exponentially large 
volume. 

recall that Gfj = 8T 8T , K and T, = n + ibi 
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Figure 1: Potential (multiplied by 10 15 ) for the canonically normalised volume modulus for a = 
0.05. 

In order to obtain a phenomenologically viable model, it is necessary to lift the AdS 
minimum to dS or Minkowski, without spoiling the stabilisation of the Kahler moduli of 
the theory. The procedure proposed in [7] is to add a D3 brane which generates a term in 
the potential of the form: 

V m *^ = e-^. (2.22) 
The full potential then takes the form 

V = (1 - a$3/2) e - V 7 ?* + ee~^. (2.23) 

By tuning e one can then find a dS or Minkowski minimum with all moduli stabilised. The 
minimum will lie in the same region as the initial AdS as illustrated in Fig. 2. 

3. Moduli/matter interaction 

The goal of this work is to investigate the interaction between matter and moduli fields. In 
particular, we want to know whether very large energy densities can destabilise the moduli 
or create bubbles of different vacua, and if so what will be the observable consequences. 

Let us first enumerate some general features of this interaction. First, the volume 
modulus couples to everything. This is primarily because it sets overall scales and enters 
all dimensionful quantities. Secondly, the volume modulus is also the lightest modulus. 
This implies that it is the most appropriate field to consider the interactions of, as it 
should be easiest to destabilise. Finally, in order to try and destabilise the field we want 
the densest possible regions in order to have maximal effect on the potential. 
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Figure 2: The full potential (multiplied by 10 15 ) for the canonically normalised volume modulus 
for a = 0.05. 



Natural candidates to consider are neutron stars. These are very dense, gravitationally 
bound systems made up mostly of neutrons. We start by noting that the energy density 
of a neutron star is: 

PNS ~ h-QCD- 

Due to gauge coupling running this term will be moduli dependent and therefore it will 
contribute to the moduli potential, Eq. (2.23). Let us briefly illustrate how this occurs. 
Recall that the QCD (3— function is 

P(g(ri) = -^9(vf, (3.2) 

where as usual /3 (<?(//)) = [i^. Substituting the definition of the j3— function into Eq. (3.2) 
and integrating one finds that 

lnW ' ) = "¥(^"^))' (3 - 3) 



and setting <?(//) — > oo, fjf = Aqcd yields 



•2 



A QC D = v"e ^ 77 >. (3.4) 
Assuming that the coupling starts to run from the string scale M s , we set 

y !' = M s = A %, (3.5) 
V v 

g 2 (v") = g 2 YM, (3-6) 



Mp 



A QCD = ^=e "Tm. (3.7) 
V V 



to find 



The crucial fact to note here is that the volume dependence comes from assuming that 
the coupling starts to run from the string scale. The physics of this is simply that the 
QCD scale is a function of the scale from which the coupling starts running. However in 
string theory with a canonical gravitational action the string scale is itself a function of 
the moduli, and so the QCD scale itself is a function of the moduli. 

This allows us to add an extra term to the moduli potential, coming from the interac- 
tion with matter 



V 

The full potential then becomes 

V = (1 - a^ 2 )e~ ^ + ee~^ + poe^, (3.9) 

where we take 

Po <xA A QCD e^*°. (3 . 10) 

At this point we must note that the assumption that the gauge coupling starts to run 
from the string scale can be relaxed. For example, in local models it is not the string scale 
but instead the winding scale from which couplings start running [11]. One can then take 

// oc V~ n , (3.11) 

which implies Aqcd oc exp(— ~ <I>), and therefore the term in the potential that 
parametrises the interaction with matter becomes 



Aqcd oc e~v» (3 . 12) 

Some potentially interesting cases include the Kaluza-Klein scale, Mkk = Mp/V 2 ^ and 
the Unification scale Mqut = Mp/V 1 ^ 3 . Throughout the rest of this paper we will study 
the case where the couplings start to run from the string scale. The results from Kaluza- 
Klein our GUT scale running will be essentially the same with the feature that the larger the 
n, the smaller the modulus vev shift and the denser the environment required to destabilise 
the modulus form its vacuum vev. 

In the interests of having a well-defined model, we will assume throughout that the 
coupling to matter has the form derived above and the potential for the modulus field is 
given by Eq. (3.9). 



4. Analysis of the volume modulus potential 

In this section we analyse the potential for the volume modulus, Eq. (3.9), investigating 
how the matter contribution can distort the potential and potentially lead to runaway. This 
will happen in the region of moduli space where the local energy density is comparable to 
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the combination of the LARGE volume potential plus uplifting term. For a fixed matter 
energy density, it is possible to achieve this by tuning the a parameter in the potential, 
which is related to the Euler number of the extradimensional manifold. 2 For fixed a, the 
same effect occurs as the matter energy density is increased. 

This potential tuning process is not completely free. One important constraint to this 
analysis comes from fifth force experiments. These limit the range of allowed masses for 
the volume modulus. For gravitational strength fifth force, the modulus mass must lie 
outside the range [KT 17 , 1(T 2 ] eV (see e.g. [12]). We will first analyse the case where the 
local source is a neutron star and then perform a more generic analysis. 

4.1 Neutron Stars 

We first examine the possibility of shifting the moduli vevs in a neutron star. We start by 
analysing the case where moduli physics is Planck coupled and then relax this assumption, 
treating the coupling as an extra free parameter. 

We can approximate the potential in the vicinity of the minimum as 

V{*) = m 2 (* - $o) 2 + A 4 QCD e-^-^/ M \ (4.1) 

Here Mj = XMp gives the coupling strength of the modulus - if A ~ then the 

modulus is Planck coupled, whereas if A <C 1 then the modulus-matter coupling is stronger 
than gravitational. <3?o is the vacuum expectation value of the modulus. 

From this we see that the shift in the modulus away from its vacuum value <l>o is given 

Recalling that in string theory the high energy couplings are directly related to the vevs of 
moduli fields, one observes that the fractional shift in a dimensionless coupling is of order 

(g-$o) _ A|gg 

M x ~ 2m? M 2 X ' y ' 

These expressions are valid for small shifts in the modulus vev: for larger shifts, the global 
structure of the potential will become relevant. 

4.1.1 Planck coupled moduli 

In general we expect moduli to be Planck-coupled, and this is true for the volume modulus 
in the LARGE volume models. Since 

A QCD ~ 10 _80 Mj,, in order to get the minimum inside 
the compact object to differ from the one outside we are required to work at very large 
where the local energy density effects become comparable to the background potential. 
This in itself does not pose any problem. However the mass of the modulus is given by Eq. 
(2.16) and is also determined by <1>. For Mx ~ Mp we see from Eq. (4.3) that an 0(1) shift 
in couplings (equivalently an 0(Mp) shift in the modulus vev) requires a modulus mass of 
m ~ 10 -11 eV (this corresponds to a ~ 3.1 x 10~ 3 and ~ 48). However this mass value 



2 This modifies the gravitino mass and gravity-mediated susy breaking scale in the theory. 
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falls within the range that is excluded by consideration of fifth force experiments and in 
the context of the large volume models would also require a string scale of around lOkeV, 
which is manifestly excluded. 

4.1.2 Strongly coupled moduli 

In a more phenomenological approach one might consider allowing the coupling strength 
Mx to be a free parameter rather than tying it to the Planck scale Mp, even though in the 
model considered the volume modulus with the specific potential, Eq. (2.23), is necessarily 
Planck coupled. Moduli in string models can certainly be coupled much more strongly than 
Mp (in the large volume models the blow-up moduli have matter couplings suppressed by 
Ms for example). 

In this case we require the modulus mass not to be smaller than 10 -3 eV and require an 
0(1) shift in a coupling. Depending on the precise value of the modulus mass, a numerical 
analysis now reveals that the interesting range for M x is between 10 7 and 10 10 GeV, 
consistent with the estimates of Eq. (4.3). Such strongly coupled moduli are in principle 
obtainable for models with low string scales. However this does not really help us, as we 
expect light moduli with masses m < leV and couplings only suppressed by Mx ~ 10 9 GeV 
to be excluded by bounds on the cooling of SN1987A by emission of exotic light particles. 3 
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Figure 3: Modulus mass as a function of 
the a parameter for, from left to right, M x = 
10 10 ,10 9 ,10 8 ,10 7 GcV 



Figure 4: Normalised A<fi as a function of 
the a parameter for, from left to right, M x = 
10 10 ,10 9 ,10 8 ,10 7 GeV 



4.2 Cosmic strings and other VERY dense objects 

We now put neutron stars to one side and consider the typical density that would be 
required in order to cause Planckian displacements of moduli. We stress that we here keep 

3 The cooling of SN1987A excludes axion decay constants f a < 10 10 GeV - see for example [13]. The 
precise numbers entering the bound do depend on the pseudoscalar nature of the axion coupling, and so 
would be modified for emission of a modulus, but we do not expect the order of magnitude of the bound 
to change. 
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the requirement that moduli physics is Planck coupled, as expected from string theory. We 
again aim to find the region of parameter space which allows for the shift of moduli vevs 
while remaining compatible with constraints from fifth-force experiments. 

The fundamental reason why neutron star energy densities could not destabilize Planck 

' ^qcd\ 2 ^ n -40 



coupled moduli is due to the hierarchy 



M P l 



10 in (4.3). It is therefore clear that 



we would need objects of higher energy density. 

In what follows we perform a numerical scan for objects of energy densities p € 
[10~ 60 , 10 _44 ]M|,;. The results of the numerical study are shown in Figs. 5 and 6. We 
require an 0(0.1) fractional displacement of the moduli from its vacuum value in a dense 
background. We plot in Fig. 5 the vacuum mass of the modulus for which this can be 
attained, for several different values of the background density. One sees that as the local 
perturbation becomes less dense, the region of parameter space where a deviation of order 
1% in the modulus vev is attained corresponds to a region where the mass for this modulus 
is smaller. Keeping in mind that the fifth force lower limit for the mass is around 10 _2 eV 
we conclude that the minimum energy density of an object capable of generating regions of 
different vacuum while still being compatible with fifth force constraints is p 
Figure 6 shows the fractional shift in the modulus vev for various energy densities. 
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Figure 5: Modulus mass as a function 
of the a parameter for, from left to right, p = 
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Figure 6: Normalised A<f> as a function 
of the a parameter for, from left to right, p = 
10" 60 , 10" 55 , 10" 50 
10" 38 M^. 
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Having performed this generic analysis we should now consider where such densities 
could come from. We first note that any such object has to be extremely dense, about 20 
orders of magnitude denser than neutron stars. For static objects a couple of possibilities 
present themselves. First, there are GUT cosmic strings. These topological defects are 
remnants of the GUT breaking which may have happened early in the history of the 
universe. Their mass per unit length p is related to the scale of the breaking of the 
symmetry that generates them a by p oc a 2 . In the case of GUT strings, a ~ 10 16 GeV 
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and therefore /i ~ 10 32 GeV 2 . For cosmological purposes the strings are taken to be one 
dimensional objects, since they have to be either closed (unstable) or infinite. This is a 
simplification and from the theoretical point of view they are known to have a finite radius 
of the order of the correlation length of the field that spontaneously breaks the GUT. 

Secondly, one could imagine dark sector analogues of neutron stars - compact bound 
objects held together by degeneracy pressure in the same way as neutron stars are, in the 
case that there existed a dark analogue of QCD with a confinement scale A^ark ^ ITeV. 
To avoid black hole formation such objects would need to be highly compact with rather 
small radii. 

However, there is no strong reason to think either of these two objects exist, although 
we will further discuss the latter case in Sec. 4.4. One case which does exist and where 
highly energy-dense regions are expected is that of singularities arising either cosmologically 
or through black hole formation. This will be analysed in the next sections. 

4.3 Cosmological singularities 

The fact that near cosmological singularities energy densities may be arbitrarily high makes 
them good candidates for the study of a coupled modulus/matter system. This means that 
not only is it possible to find destabilise a minimum for the moduli potential but also that 
decompactification might be possible. 

The model we consider is that of a closed, matter dominated FRW universe. The 
overall spacetime structure is TZ + x S 3 x CY 6 . As before we assume that the matter is 
baryon-like, with a mass depending on the volume modulus as in Eq. (3.7). We take 
the modulus scalar field potential to be given as before by the LARGE volume potential 
plus the up lifting term. We assume the volume modulus is initially at the minimum and 
unexcited, so that Q^ t i n u = (p<s>/ptot)init = while £l ma tter,init = 1- The matter-moduli 
coupling comes from the field dependent mass in Eq. (3.7), and takes the form 

p(a,*) = ^e-^(*-*°>. (4.4) 
a 6 

The evolution of the system is determined, as usual, by the Friedmann equation 

(i)' = i(T + "W + '(-«)"^ <4 ' 5) 

and the Klein- Gordon equation for a homogeneous and isotropic field 

$ = 3-6 + (Wa,$) + p(a,$)) $. (4.6) 
a 

To study the coupled modulus-matter dynamics we solve the equations of motion 
numerically using a modified version of [14]. As expected the initial evolution of this 
universe is the same as a regular closed matter dominated universe: the scale factor grows 
to a maximum before the universe will start collapsing, with the scale factor shrinking to 
zero. This behaviour is shown in Fig. 8. During most of the collapse, the system will not 
show any departure from the behaviour of a closed, matter dominated FRW universe. It is 
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only when the scale factor becomes sufficiently small (let's call it a), that the term in Eq. 
(4.4) will become large enough to modify the structure of the modulus potential, Eq. (3.9), 
and play a significant role in the evolution of the field. This happens, to first order, when 
the matter contribution to the potential at the minimum is comparable to the potential 
plus matter at the maximum, i.e. 

mino maxn ■ 

(4.7) 

Note that mino an d maxQ denote the minimum and maximum of the LARGE volume plus 
uplifting potential (which differ from the extrema of the full potential). One can get an 
semi analytical estimate for a, by solving Eq. (4.7), finding 

\ 1/3 



V{^max ) 

This estimate can be refined by expanding a c = a + 5a and solving the condition 

V($) +p{a = 6a,$)\ min = V(<5>) + p{a + 8a,$)\ max , (4.9) 

where min and max denote the position of the minimum and the maximum of the full 
potential, i.e. including matter contribution, when a = a. Solving Eq.(4.9) one finds 

/ \ 1/3 

a c = a + 5a=( — ( e -V6* roi „ _ e -Ve*m«) ] ; ( 410 ) 



max) Vi'&min) 

which is in good agreement with numerical estimates. 

As the scale factor approaches a c , the modulus expectation value begins to shift and 
the barrier to decompactification decreases. After a short period of evolution, the term 
(4.4) becomes the dominant term in the potential and there is no obstacle to prevent the 
field from rolling towards infinity. The evolution of $ is shown in Fig. 8 and we see the 
sharp increase in <3? beyond a critical time. This final stage of the evolution sees a rapid 
runaway of to infinity while the scale factor continues to shrink to zero. The energy 
density of the universe will quickly be dominated by the kinetic energy of the volume 
modulus, with the dust and potential contributions becoming subdominant. The evolution 
of the system is depicted in Figs. 7, 8 and 9. 

In this last region it is possible to solve the equations of motion analytically. The 
equations of motion are 



-H = -- + ~ 2 - (4-12) 



dt\ J \d$ 

Jt ~ "T + ^ 

Once starts to runaway, the potential vanishes exponentially fast and ceases to be a 
significant contribution. In this asymptotic regime we can then solve these equations by 

Hi) = $o - y|ln (to - t) , 

/o\ 1/6 

a(t) = (2) ^- t ) 1/3 - ( 413 ) 
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Figure 7: Time evolution of the scale factor Figure 8: Time evolution of the modulus in 
in a closed, matter dominated FRW universe. a closed, matter dominated FRW universe. 
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Figure 9: Late time evolution of the energy densities of matter and scalar field (kinetic and 
potential). Dashed line tldust, dashed-dotted line Q po tentiai, full line ftkinetic- 

From the solution we see that we can self-consistently neglect the curvature term in Eq. 
(4.12) for ^± <^ l. We can relate the canonically normalised field <3?(t) to the volume of 
the compact space using Eq. (2.20), to find that the compact volume V evolves as 

V = -^- (4.14) 
to — t 

and so diverges as t — > to. 

In the limit as to - * t this therefore describes a universe where 



1. The 4-dimensional scale factor shrinks to zero size as (to — t) 1 / 3 . 



2. The volume of the six compact dimensions is first destabilised before diverging as 
(to-*)" 1 - 

3. There is a singularity at t = to, at which the 4-dimensional scale factor is formally 
zero and the 6-dimensional volume is formally infinite. 

This represents a universe which bounces, albeit heterogeneously: as the three spatial 
dimensions collapse, the extra six dimensions expand and reach infinite volume in finite 
time. 

Let us consider the validity of a 4-dimensional effective field theory treatment. From 
the solution for in Eq. (4.13) it is clear that the (kinetic) energy density of the $ field 
diverges in finite time, 

* = \* 2 = vth?- (4 - 15) 

Furthermore as decompactification occurs and the extra-dimensional volume increases, the 
4-dimensional string scale decreases. 

Mt M =^ = Ml>{t Q -tfe-^. (4.16) 

Consequently, independent of initial conditions, as t — > to the system will evolve to a state 
where the kinetic energy density in the field $ is greater than the apparent cutoff of a 4d 
effective field theory. 

Fortunately it is easy to understand what is happening from a higher dimensional 
perspective. Once the field $ gets over its decompactification barrier, its potential soon 
becomes negligible - we can see in Fig. 9 how the kinetic energy of $ is the dominant 
contributor to the energy density. Although we have written the problem in the language 
of 4-dimensional effective field theory, <I> is originally the volume modulus of the extra 
dimensions. Neglecting the potential energy, the system is then fundamentally that of 10- 
dimensional general relativity. The above solution then corresponds to a Kasner solution 
dimensionally reduced to 4 dimensions. 

It is not immediately clear that the above numbers are consistent with a Kasner solu- 
tion. Recall that the Kasner solution is 

ds 2 = -dt 2 + ^t 2p >dx 2 , (4.17) 

i 

with ^2 Pi = 1 an d 2 Pi = 1 • For a 1 + 3 + 6 dimensional Kasner solution, the allowed 
exponents are (using p% to denote the 3-dimensional growth and p$ for the six-dimensional 
growth) 

11 5 —1 

P3 = -g,P6 = g, P3 = g,P6 = -^-- (4.18) 

It is clear that the evolution of the scale factor in (4.13) and the volume in (4.14) do not 
fit these conditions. However, note that the metrics are different: the metric used in the 
Kasner solution (4.17) is the 10-dimensional string frame metric, whereas the 4-dimensional 
metric for which Eqs. (4.13) and (4.14) applies is a dimensionally reduced metric that is 
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related to the 10-dimensional metric by factors of the internal volume. It is then expected 
that the 4-dimensional scale factor does not have the Kasner exponent appropriate for a 
10-dimensional (1+3+6) solution. 

There is one striking feature about this behaviour. The initial conditions (a closed 
matter-dominated FRW universe) were unexceptional. However these conditions unavoid- 
ably evolve to give dynamic super- inflationary behaviour of the compact dimensions, which 
in theory reach infinite volume in finite time. Furthermore, this behaviour commences in 
the region controlled by effective field theory, where a 4d description is valid. In prac- 
tice, the evolution of a Kasner solution should break down as the contracting dimensions 
approach the lOd string scale, regulating the infinity. This super-inflationary behaviour 
cannot be prevented - growth in the 4d energy density is a necessary consequence of a 
spacetime crunch - and this energy density must always eventually overcome the barrier to 
decompactification. 

As formulated the dynamics have started with 3 large and 6 compact dimensions, 
ending with 6 large and 3 compact dimensions. However the physics is such that there 
is no reason not to reverse the process, and imagine starting with 6 large and 3 compact 
dimensions and ending with 6 compact and 3 large dimensions. In effect the universe 
bounces within effective field theory, which is achieved by the bounce occuring in differ- 
ent dimensions to the collapse: the collapse of certain dimensions triggers the expansion of 
others. This physics has some similarities to pre-Big Bang cosmology [15, 16, 17]. This was 
formulated using the heterotic dilaton and the O (d, d) symmetries of toroidal compactifi- 
cation. It would be interesting to make the connections more precise and see whether this 
super-inflationary growth of the compact dimensions is able to mimic some of the physics 
of conventional inflation. We leave this to future work. 

4.4 A static solution 

In this section we study spherically symmetric configurations and analyse the resulting 
modulus profile. The models considered in this section are similar to the 'dark stars' 
mentioned in Sec. 4.2: we look for stable solutions of matter coupled to a modulus field, 
protected against gravitational collapse. In Sec. 4.5 below we allow for dynamical evolution 
of the profile and study gravitational collapse of dust balls. 

Consider a spherical ball of dust of radius R. The dust has energy density p but is 
pressureless (P = 0). Let the line element inside the dust ball be curved FRW: 



where a is the scale factor and k is the spacial curvature. By Birkoff's theorem, the 
spacetime for r > R is the Schwarzschild solution: 




(4.19) 




r 



(4.20) 



where M is the total mass in the region r < R. 
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For arbitrary values of p and P, the equations of motion for the scale factor a(t) in 
Eq. (4.19) are: 



° = -(p + 3p), (4.21) 
a 

, 2 

a \ p k 
= 3 ~ o2' 



(4.22) 



One can find a static solution (d = a = 0) of Eqs. (4.21), (4.22) by taking 



k = p/3, (4.23) 
P = -3p, (4.24) 

where we have set a = 1. This is the Einstein Static Universe. 

In the spirit of the previous sections, let the mass of the dust particles be a function 
of the modulus field <j). The dust energy density is therefore given by: 

Pdust = Poe'^t (4.25) 

which represents a source for the volume modulus potential in the region r < R. The 
dust has zero pressure but there is a pressure contribution from the vacuum energy of the 
modulus field as it is displaced from its minimum. 

The simplest modulus profile compatible with the equations of motion is obtained by 
considering three distinct regions. The first region is r < R, in which the field is at the 
minimum of the effective potential V e ff = V& + pdust- For r > R + l/rru the field is at 
its vacuum minimum. In the transition region R < r < R + 1 jm^ the field interpolates 
smoothly between the two distinct minima. 

The volume modulus profile described above implies that inside the dust ball, 

P = Pdust + Vj,, (4.26) 
P = -V h (4.27) 

where these are evaluated at the minimum of the combined potential V^ust + V<j>- Although 
there is no pressure from the dust, the displacement of scalar fields from the minimum leads 
to a contribution to vacuum energy. One then finds that Eqs. (4.23) and (4.24) become 

k= Pdust + (428) 

3 

Pdust - 2V^ = 0. (4.29) 

Since one can treat the spatial curvature as a free parameter, the condition for existence 
of a static solution for the region r < R reduces to Eq. (4.29). As long as the combined 
potential V = Vdust + ^> exhibits a minimum where Vdust = 2V^, then a static solution 
will exist. There are two possibilities of tuning the system to generate a solution to Eq. 
(4.29): one may either tune the dust density po or the parameter a in the LARGE volume 
potential. There is then effectively a 1-parameter set of solutions parametrised by the 
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density of the interior. There are various constraints on this set, for example by limits on 
the volume modulus mass coming from fifth force experiments, as discussed in Sec. 4. 

In Table 1 we display four different static solutions. These were obtained by fixing 
a and then tuning the density p. It is interesting to study the properties of the dust 
distribution that sources the nontrivial volume modulus profile, in particular its mass and 
radius. Given that the density is uniform and fixed by Eq. (4.29), the mass will be given 
by M = ^pR 3 , where R is the radius of the spherical dust distribution. One may write 
the radius in terms of the Schwarzschild radius as R = £,R sc h, where £ > 1. Since, by 
definition, R sc h — 2GM one finds 



R 



this implies that 



sch 



R 



lirGpt 3 



(4.30) 



(4.31) 



8irGp£ 

One then concludes that the properties of the dust ball are completely determined by 
£ since p is fixed by requiring a static solution. Note that the smaller £ the larger the star's 
radius and mass. In Table 1 we display the radii and masses for the four cases under study 
considering £ = 1.1. 



a 




P«) 


RiM- 1 ) 


M(M p ) 


0.05 


1.42 x 10' 


-7 


5.39278 x 10' 


-7 


449 


204 


0.01 


4.45 x 10- 


-19 


5.36521 x 10- 


15 


4.49 x 10 6 


2 x 10 6 


0.005 


2.80 x 10- 


-29 


6.288865 x 10 


-22 


1.3 x 10 10 


6 x 10 9 


0.001 


5.16 x 10- 


-82 


2.15993 x 10- 


57 


7 x 10 27 


3 x 10 27 



Table 1: Mass of the canonically normalised volume modulus, density, radius and total mass of 
the dust sphere, as functions of the a parameter in the LARGE volume potential. Radius and total 
mass computed assuming £ = 1.1. 

The numerical results in table 1 reveal that there is a very large hierarchy both in 
radius and mass between the various cases studied here. This exemplifies the issue raised in 
section 4.2: for masses of <3? greater than that allowed by fifth force constrains, namely m$ > 
10- 30 M P , the size of such objects is extremely small (R < 10 23 cm). The corresponding 
mass is M < 1kg. While they could in principle form part of dark matter, it is hard 
to see how interesting physics can be extracted from them. The local density of dark 
matter objects with kg masses is not larger than 10~ 21 m~ 3 and so objects would be both 
unobservable and undetectable. 

4.5 A dynamic solution: black hole formation 

In this section we generalise the analysis of Sec. 4. 4 to allow for dynamical evolution of 
the system. While Sec. 4.4 was restricted to static solutions, here we look for collapsing 
solutions of the coupled matter-modulus system. If one lets the gravitational collapse 
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last for long enough, the final state will be a Schwarzschild black hole. However before one 
reaches the singularity, the local density becomes arbitrarily large, which could be sufficient 
to destabilise the modulus. 

Here our initial conditions are a large, dilute dust ball which we allow to collapse 
towards a black hole. As before we assume that the spacetime is given by a curved FRW 
universe smoothly connected at the surface of the dust sphere (r = R) to a Schwarzschild 
solution. In contrast to Sec. 4.4, we now use the Friedmann equation to obtain dynamical 
solutions of the scale factor and of the energy density for r < R. 

The profile for the volume modulus is obtained by considering: 

f0 + 3|0 + V^ = O ,re[0,R] 
<f>=l (1 - 2M/r)<j>" + (f(l + 2M/r) + 2M/r 2 )^ - = ,r € [R,R + m~ l ] (4.32) 
[4>oo ,r € [R + m^ 1 ,oo] 

where ' = , ' = , M is the mass sourcing the Schwarzschild geometry and <poo is the 
value of the field in the minimum of its vacuum potential. 

The time dependence of the system arises in two different ways. Inside the dust 
distribution one has a homogeneous and isotropic positively curved universe, which is in 
general dynamic. This causes the scale factor and the volume modulus to be functions of 
time. For r > R we assume that the time dependence comes only from the time variation 
of the radius of the dust ball (which can be traced to the time variation of the scale factor). 
In particular, we assume that the scalar field profile adjusts instantaneously to the change 
in dust density, i.e. that the characteristic time scale of the gravitational collapse is much 
larger than the typical timescale of the variation of the scalar field. 

One must clarify what is meant by M in Eq. (4.32). Naively one would expect the mass 
of the dust ball to be given by the volume integral of the dust energy density. However there 
are also non negligible contributions from the potential and kinetic energies of the volume 
modulus. These are, for most of the evolution of the system, one order of magnitude smaller 
than the dust contribution but become more significant as the system evolves. Taking this 
effect into consideration, the mass sourcing the Schwarzschild geometry is 

M = J Sxyfg (p dust + p vw + p K (4,)) ■ (4.33) 

In Fig. 10 we compare Eq. (4.33) with the naive estimate. The deviation between the two 
estimates is negligible throughout most of the evolution but increases with time, diverging 
at the end. This divergence is due to the fact that when the density reaches a critical 
value, the volume modulus potential no longer exhibits a minimum and decompactification 
happens. This means that the modulus starts to roll and its kinetic energy dominates the 
energy density inside the dust ball, in a similar way to the behaviour described in Sec. 4.3 
for cosmological singularities. 

The time evolution of the system is depicted in Fig. 11. We start with a spherical 
distribution of radius r > 2M, regime I. In this regime the whole transition region (shaded 
area) is outside the horizon and the deviation of the volume modulus from its vacuum 
minimum could in principle be observable. As the gravitational collapse evolves a black 
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Figure 10: Comparison between the naive estimate for 2M (dashed line) and the result of Eq. 
(4.33) (solid line). The shaded area represents the region where the field varies between the value 
inside the dust ball and the vacuum minimum. 

hole will form and subsequently the transition region will start to fall inside the horizon, 
this is regime II. In this regime, it is still in principle possible to observe the consequences of 
the matter modulus interaction since part of the transition region lies outside the horizon. 
When the whole of the transition region is veiled by the horizon, regime III, an observer 
sitting outside the black hole will not be able to measure the nontrivial profile of the volume 
modulus. 




Figure 11: Sketch of the black hole formation process. The dashed line depicts the Schwarzschild 
horizon and the shaded area depicts the region R < r < R+ 1/ where the field rolls between its 
value inside and outside the dust ball. 

The time evolution of the volume modulus profile in regime I is shown in Fig. 12. As 
the system collapses and the density increases, the volume modulus vev inside the dust ball 
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Figure 12: Time evolution of the volume modulus profile in regime I. 



increases and the transition region where the field is allowed to vary moves to the smaller 
radius region. The evolution of the field profile in regimes II and III is qualitatively similar 
to the one just described. 4 

5. Discussion 

The purpose of this paper has been to consider the interaction of moduli and matter fields, 
and specifically to analyse the circumstances under which moduli fields can be destabilised 
from their vevs by dense concentrations of matter. We have described the origin of mod- 
uli/matter couplings and the assumed form of the moduli potential. We have tried to 
consider 'honest' values for parameters such as coupling strengths and moduli masses, 
enforcing the consistency constraints that emerge from string compactifications. 

Our results have both negative and positive elements. On the negative side, it is not 
feasible to destabilise moduli through even the densest astrophysical environments. We 
considered neutron stars as the densest known form of matter, and found that Planck 
coupled moduli could only be destabilised if their masses were deep in the region excluded 
by fifth force constraints. If the coupling is relaxed from Planck strength, then the values 
of the coupling for which destabilisation is possible are excluded by the stellar cooling 
constraints that require f a < 10 9 GeV. This only leaves very exotic cases such as cosmic 
string and hypothetical dark analogues of neutron stars. 

On the positive side, our results show that the modulus/matter coupling will play a 
significant role in gravitational collapse. The cosmological collapse of a matter-dominated 
FRW universe was shown to lead to a super-inflationary decompactification of the internal 

4 In regime II one must impose fmiteness and continuity of the solution at the horizon 
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dimensions as the volume modulus is destabilised. This super-inflationary decompactifi- 
cation leads to the internal dimensions reaching infinite volume at the same time as the 
external scale factor vanishes. It would be interesting to see whether this super-inflationary 
decompactification can mimic some of the features of conventional inflation and we leave 
this for future work. In a similar vein we also saw that the process of black hole formation 
will lead to decompactification during the period of collapse. 
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